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Abstract 


We prove an effective integrability criterion for differential-algebraic 
Pfaffian systems leading to a decision method of consistency with a 
triple exponential complexity bound. As a byproduct, we obtain an 
upper bound for the order of differentiations in the differential Null- 
stellensatz for these systems. 
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1 Introduction 

Let x := xi, ..., x m and y := y\,... ,y n be two sets of variables; the first 
ones represent the independent variables (i.e. those defining the partial 
derivations) and the second ones are considered as differential unknowns. 

The notion of Pfaffian system is introduced by J.F. Pfaff in [20] and in 
its simplest form it is defined as a system of partial differential equations of 
the type: 



(1) 
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where each fij is an analytic function around a certain fixed point (xo, yo) € 
C m x C n . 

The properties of these systems were extensively studied during the 
XIXth century by notable mathematicians as Jacobi [11], Clebsch [2] and 
Frobenius [ 6 ] (see also [10] for a detailed historical approach). One of the 
main problems considered is the so-called complete integrability of a Pfaffian 
system X: the existence of a neighborhood U of the point (xo, yo) such that 
for all point (x, y) £ U there exists a solution 7 of X such that y(x) = y. A 
complete solution of this problem is known today as the “Frobenius Theo¬ 
rem” (see [ 6 ]): 

Frobenius Theorem Let X be the Pfaffian system (1). Then X is com¬ 
pletely integrable in (xo,yo) if and only if for all indices i,j,k with i = 
1 ,n, and j,k = 1 ,... , m, the function Dj(fik) — Dk(fij) vanishes in a 
neighborhood 0 /(xo,yo)- 


Here, Dj (h) denotes the j-th total derivative with respect to X of any an- 

j y — f ■ 
^ 1 dyi * 13 ' 

Observe that the vanishing of the functions Dj(fik) — Dk(fij) is clearly a 
necessary condition because of the equality of the mixed derivatives of an 
analytic function, but the sufficiency is not obvious. 


dh 

alytic function h in the variables (x, y), namely Dj(h ) := ——- 


We remark that Frobenius Theorem is more general than the statement 
above because it remains true also for systems of differential linear 1 -forms 
(today called Pfaffian forms). In this sense Frobenius’s article may be con¬ 
sidered as a main source of inspiration for the transcendental work by Cartan 
about integrability of exterior differential systems developed in the first half 
of the XXth century (see [1]). 

In Cartan’s theory two notions play a main role: the prolongation and 
the involutivity of an exterior differential system. In the particular case 
of differential equation systems these notions can be easily paraphrased. 
The prolongation of a system consists simply in applying derivations to 
it. On the other hand, an involutive system is a system with no hidden 
integrability conditions; in other words, no prolongation is necessary to find 
new constraints. Roughly speaking, the involutive systems are those to 
which the classical method of resolution by means of a power series with 
indeterminate coefficients (known as Frobenius method) can be applied in 
order to search for a solution or to decide that the system is not integrable. 

A main general result related to the integrability in Cartan’s theory is 
the Cartan-Kuranishi Principle: Any (generic) exterior differential system 
can be reduced to an equivalent involutive system by a finite number of pro¬ 
longations and projections. This result, conjectured by Cartan, is finally 
proved by Kuranishi in [16] (see also [18]). Despite the algebraic precisions 
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given during the ’60s and ’70s (see for instance [24, 9, 23]) and more recently 
in [ 21 ] and [ 22 ], up to our knowledge, no effective a priori upper bounds for 
the number of required differentiations nor projections are established. 

In this paper we use arguments close to the Cartan-Kuranishi Principle 
(namely, prolongations and projections) in order to study the integrability 
(not necessarily complete) of differential-algebraic Pfaffian systems. More 
precisely, a differential-algebraic Pfaffian system is a system of partial dif¬ 
ferential equations X as follows: 

E:=(i =f<X - y) ( 2 ) 

l g(x,y) = 0 

where f := ( fij)ij and g := g\,... ,g s are polynomials in C[x,y]. As before, 
x := (xi,..., x m ) and y := (yi, ..., y n ) denote sets of variables. 

By means of prolongations and projections we construct a decreasing 
chain of algebraic varieties C m+n 2 Vo 3 Vi 3 which becomes station¬ 
ary at most at the (m + n + l)th step. If we denote by Vqq the smallest 
variety of the chain, we prove the following criterion for integrability (see 
Theorem 15 and the comments that follow it): 

Theorem 1 The system X is integrable if and only if the algebraic variety 
Voo is non-empty. The analytic variety of all the regular points ofY 00 is the 
biggest analytic variety containing all the graphs of analytic solutions o/X. 

Moreover, by means of basic algorithms from commutative algebra, the 
criterion can be transformed in a decision algorithm which runs within a 

complexity of order ( nm.ad) 2 ° {n+m) , where d is an upper bound for the de¬ 
grees of the involved polynomials and a := max{l,s} (see Theorem 19 ). 
As a byproduct, we obtain an upper bound of the same order for the order 
of differentiations in the differential Nullstellensatz for differential-algebraic 
Pfaffian systems (see Theorem 23). In this sense, the present work can be 
seen as a continuation of [3]. 

The paper is organized as follows. In Section 2, first, we introduce the 
notation we use throughout the paper; then, we show how to reduce the inte¬ 
grability problem of general (non-autonomous) differential Pfaffian systems 
to the autonomous case; finally, we prove our integrability criteria. Section 
3 is devoted to analyzing quantitative aspects of the problem: we present an 
effective decision method for the integrability of differential-algebraic Pfaf¬ 
fian systems and we prove an upper bound for the order in the differential 
Nullstellensatz for these systems. 
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2 A geometrical criterion for the integrability of 
differential-algebraic PfafRan systems 

In this section we exhibit a necessary and sufficient criterion for the inte¬ 
grability of differential-algebraic Pfaffian systems in terms of the dimension 
of a finite decreasing sequence of algebraic varieties associated to the dif¬ 
ferential system (see Definition 7 and Theorem 15). These varieties can 
be constructed explicitly by means of three basic operations: prolongations 
(i.e. differentiations), linear projections and reductions (i.e. computation of 
radicals of polynomial ideals). 


2.1 Notations 


Let m, n G N. We consider two sets of variables x := aq,..., x rn (the so- 
called independent variables ) and y := yi, ..., y n (the variables playing the 
role of differential unknowns ). For each pair 1 < i < n, 1 < j < m, let 

fij be a polynomial in the ring C[x, y]. Finally, let g := g \,..., g s € C[x, y] 
be another finite set of polynomials. 

A differential-algebraic Pfaffian system is defined as a partial differential 
system £ of the type: 


£ = < dx __ 

5fc(x, y) = o 

or in its simplified form: 


dv' 

■KT = /b'( x >y)> 1 < i < n, 1 < j < rn, 


1 < k < s, 


£={ ll =f<x ' y) ’ 

l g(x,y) = 0 


( 3 ) 


where f denotes the set of polynomials fij. If the polynomials f and g do not 
depend on the variables x, we say that the Pfaffian system £ is autonomous. 

For each index j = 1,,m and any polynomial h G C[x, y] we define 
the total derivative with respect to xj induced by £ as 


Dffh) := 


dh ^ dh 
dx ^ 




Observe that Dj(h) belongs to the polynomial ring C[x, y] and that the 
operator Dj is a derivation in this ring. 


2.2 Frobenius compatibility conditions 

We are interested in the integrability (or solvability ) of differential-algebraic 
Pfaffian systems. This notion should be understood as the existence of an 
analytic solution 7 defined in an open subset of C n with target space C m . 
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A stronger notion is that of complete integrability , which means that 
there exists a neighborhood U C C m+n around (xo,yo), such that for any 
(x, y) € U there exists a solution 7 of X verifying y(x) = y (note that this 
notion only makes sense in the case that no algebraic constraints appear, 
because such a constraint never contains a nonempty open set). If this is 
the case for the system X, the classical Frobenius Theorem (see for instance 
[4, Ch.X, §9]) gives a simple criterion for complete integrability: X is com¬ 
pletely integrable around (xo,yo) if and only if the compatibility conditions 
Dj(fie) — De(fij) = 0 hold for all triplets i,j,£ in a suitable neighborhood of 
(x 0 ,yo)- 

Our systems are in a certain sense more general than those considered 
in Frobenius’s result (more constraints appear) and we are interested in a 
weaker notion of integrability. In order to begin our analysis, we introduce 
the following ideal which will enable us to deal with the Frobenius compat¬ 
ibility conditions and will be involved in the first step of our criterion: 

Notation 2 For a given Pfaffian system X as in (3) we denote by $ the ideal 
in the polynomial ring C[x, y] generated by the polynomials Dj{f i k) — D} i (fij ) 
for all indices i,j,k. 

If no algebraic constraints appear in X, the classical Frobenius Theorem 
in the differential-algebraic setting can be restated as follows: X is com¬ 
pletely integrable if and only if the polynomial ideal $ is zero. On the other 
hand, if such a system X is integrable, the ideal ^ must be properly con¬ 
tained in C[x, y], since if 7 is a solution, all polynomial in J vanish at the 
vector y(x). However, this condition (namely, the properness of 5) is not 
enough to guarantee simple integrability: 


Example: Consider the Pfaffian system 


dy = . 2 9 y_ 

dx\ ^ ’ 8x2 


2/ 2 Al¬ 


in this case, n = 1, m = 2, and no constraint appears. It is obvious that if 
7 is a solution, then it is not a constant function because the polynomials 
y 2 and y 2 + 1 have no common zeros. The ideal $ is generated by the 
polynomial 2 y(y 2 + 1) — 2 yy 2 = 2 y and then it is proper. On the other 
hand, any solution 7 must verify the equation 2y = 0 and in particular 7 is 
a constant, leading to a contradiction. 
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2.3 From non-autonomous to autonomous systems 

We begin by showing that the problem of the integrability of a (general) 
differential Pfaffian system 


£ = 


§ = f(x ' y) 
g(x,y) = 0 


can be reduced to analyzing the integrability of an autonomous Pfaffian 
system. 

To this end, we transform £ into an autonomous Pfaffian system £ aut in 
the obvious way, by introducing new differential unknowns w := w\,..., w m 
(as many as independent variables): 


£ant — < 


dy 

<9x 


= f(w,y) 


<9w 

g( W ,y) = 0 


( 4 ) 


where 5 denotes the Dirac-symbol (i.e. Sjg = 1 if j = £, and 0 otherwise). 
Clearly £ aut is an autonomous system in m+n unknowns and m independent 
variables; moreover, it is equivalent to £ from the integrability point of view: 


Proposition 3 The Pfaffian system £ is integrable if and only if the au¬ 
tonomous Pfaffian system Ti au t is integrable. 

Proof. Let 7 : U C C m —>• C n be a solution of £; then if : U —>• C m+n 
defined as V’( x ) := (x, 7 (x)) is a solution of £ aut . Reciprocally, let if : 

V —>• C m+n be a solution of £ aut in a connected neighborhood of a certain 
xi € C m . Prom the second part of the defining equations of £ au t one 
deduces that if(x.) = (x + A,^(x)) for a suitable vector A € C m and an 
analytic function £ : V —>• C n . Take xo := xi + A and IA := V + A. Define 
7 (x) := f(x — A). Thus 7 : U —> C n verifies: 

= ° lr " l+ o* —— = /b(^( x “ A )) = /ij(( x -A)+A,^( x -A)) = f i:j (x, 7(x)) 

and 

g( x > 7( x )) = g(( x - A) + A, f (x - A)) = g (if(x - A)) = 0. 

In other words, 7 is a solution of £. ■ 
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2.4 Prolongation chains 


Taking into account the result of the previous section, for simplicity, we 
restrict now our attention to considering the autonomous case, namely, 
differential-algebraic Pfafhan systems of the form 


E = 


g(y) = o 


( 5 ) 


where f and g are finite sets of polynomials in C[y]. In this setting, for a 
polynomial h G C[y], the total derivative of h with respect to xj induced by 
E is 

71 r)h 

eCly], 

i =1 

First, we introduce the so-called prolongation of an ideal in C[y]: 


Definition 4 Let I C C[y] be an arbitrary ideal and F be a finite system 
of generators of I. We denote by I the ideal of the polynomial ring C[y] 
generated by I and the polynomials Dj(f) for all f G F and j = 1,..., m. 

Proposition 5 The ideal I is independent of the system of generators F of 

I. 


Proof. Let G be another system of generators of I. It suffices to see that 
for any k = 1 ,,m and any g G G the polynomial Dffg) belongs to the 
ideal generated by I and the polynomials Dj(f). Writing g = E Pff and 


using that Dj is a derivation in C[y], we deduce: 


/e F 


d m = Y Djipf n = D j(pf) ■ f + Y p r D otf)- 

/6F /<= F /GF 


Since the first term belongs to (F) = I and the second one to the ideal 
generated by the polynomials Dj(f), the proposition follows. ■ 


We start proving an elementary test to check integrability by means of 
a prolongation and a projection, which will be the key result to our criteria. 

Lemma 6 Let E be an autonomous differential-algebraic Pfaffian system: 

£ = ( i = f(y) - 

l g(y) = o 

where f and g are finite sets of polynomials in C[y]. Let T := sj (g) and 

J := vj + Z, where 5 C C[y] is the ideal generated by the Frobenius con¬ 
ditions, as in Notation 2, and let V D W be the varieties defined by T and 
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J respectively. Suppose that V and W are the same nonempty variety in 
a neighborhood of a common regular point Q € C n . Then £ is integrable. 
Moreover, there exists a solution of £ passing through the point Q and con¬ 
tained in W. 


Proof. From the hypotheses, the ideals I and J define the same nonempty 
algebraic variety C locally around Q. Let go be a finite system of generators 
of X. Since Q is a regular point of C, we can apply the Theorem of Implicit 
Functions to the system go = 0 and so, without loss of generality (reordering 
the variables if necessary) we may suppose that there exists an integer r > 0 
such that, if y := 2 / 1 ,..., y r and y := y r +i, ■ ■ ■, y n , the variety C is locally 
defined around Q = (y 0 , yo) as the graph of an infinite differentiable function 
c p defined in a neighborhood U C C r of y 0 with target space C n ~ r . The 

£ go 

integer r is the dimension of C, the rank of the Jacobian matrix —— in a 

dy 

neighborhood of Q in C is constant and equal to n — r, and moreover, the 

Qgo 

submatrix —— has locally maximal rank n — r. 

o y 

Consider the following partial differential (not necessarily algebraic) Pfaffian 
system £ induced by £ and (p: 


£ : = 


dy 

<9x 


= f(y,<p(y)) 


where f := ( fij) with i = 1 ,..., r, j = 1 ,..., m. 

The proof of the lemma will be achieved by showing the following facts 
concerning the system £: 


• Claim 1. The system £ is completely integrable in a neighborhood of 

yo- 

• Claim 2. If y is a solution of £ around yo, then 7 := (/r, pop) is a 
solution o/£. 


Proof of Claim 1. From Frobenius Theorem it suffices to prove that for each 
pair j, k, l < j, k < m, and alH = 1 ,..., r, the relation 



dfiji y,p(y)) 

dyh 


fhk 


dfiki y,y(y)) 
dyh 



= 0 

(y,¥>(y)) 


holds for all y € U (where U a suitable neighborhood of y 0 ). 

From the Chain Rule, for any q = 1,..., m, and i,h = 1,..., r, we have the 
equality 


dfi q {y,p(y)) 

dyh 


dfi 


iq 


dyh 


+ £ 

(y>v(y)) i = r +1 


dfiq 

dye 


— (y) 

(yMy)) 
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where ip r + 1 , • • • ,f n are the coordinates of the function ip. 

These equalities can be collected in a unique relation in terms of a product of 
matrices evaluated in (y, ifi(y)) which describes the Frobenius compatibility 
conditions for the system £ as follows: 


dfij dfij_ _ £hp\ 
dy dy dy) 


7 fik ^ 

( dfik 

dfik 

dip\ 

( hi 

l frk j 

V dy 

dy 

dy ) 

\ hj 


( 6 ) 


dip 

where 7 — denotes the 
dy 

to the variables y. 


(n-r)xr Jacobian matrix of the map ip with respect 


On the other hand, differentiating the identity go(y,ip{y)) = 0 and multi¬ 
plying by the column vector we infer that, for all j = 1 ,..., m, 


dgo 

dy 


( h3 \ 

V frj j 


dgo dip 
dy ' dy 


( hi 

\ frj 


= 0 


(7) 


holds at all the points (y, ip(y)). 

Finally, from the hypothesis C coincides with V and W in a neighborhood 
of the point Q, we may suppose (shrinking the neighborhood U of y 0 if 
necessary) that h(y, tp(y)) = 0 for all h € J, in particular since Dj(g ) € J 
for all j = 1 and all g G go we have 


dgo 

dy 


hi 

V fnj 


0 


in any point (y, <^(y)), for all index j. Again, separating the variables y as 
before, this equality gives 


dgo 

dy 


hi 

\ frj 


dgo 

dy ' 


7 f(r+l)j 

V fnj 


0 . 


( 8 ) 


Q go 

Comparing equalities (7) and ( 8 ) and recalling that the matrix has full 
rank n — r we infer that 


dip 

df 


hi 

\ frj 


7 f(r+l)j 

V fnj 


in any point (y, <^(y)) and for all index j. 


(9) 
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In order to finish the proof of Claim 1 it suffices to prove the validity of 
relation ( 6 ). Since, by definition, the Frobenius conditions associated to the 
system X are contained in the ideal J, we have that 


9/ij ( / “ 

% ' l /l* 


Bfik ( Al ' 

ay ' l L J 


holds in a neighborhood of Q relative to the variety C. Splitting the variables 
y into y and y as before, this relation can be written 


dfij 

dy 


( flk 
V frk 


dfij 

dy 


f(r+l)k 


fnk 


dfik 

dy 


( h3 \ 

\ frj j 


dfik 

dy 


^ f(r+l)j 
\ fnj 


Identity ( 6 ) is obtained from this formula simply by replacing with relation 
(9) for the indices j and k. This finishes the proof of Claim 1. 

Proof of Claim 2. Suppose that p is a solution of X around y 0 . Clearly the 
image of 7 := (p, (p o /f) is contained in C because it is included in the the 
graph of ip. In particular 7 verifies the algebraic constraint g = 0. 
Decompose 7 = ( 7 , 7 ) where 7 := p and 7 := ip o p. For all j = 1,..., m we 
have the relations 


and 

dy 

dxj 


dy _ dp 
dxj dxj 


( k) \ 

\ f rj (ji,(p op) / 


( hj{l) 

\ frj( 7 ) 


d(p°p) 

dxj 


dp dp 
dy dxj 


dip 

dy 


( fij(k,<P°k) 

\ frj(p,p°p ) 


^ /(r+l)j'(T) 

\ fnj{ 7 ) 


where the last equality follows from identity (9). This completes the proof 
of Claim 2 and therefore also of the lemma. ■ 


2.5 An integrability criterion 

We introduce an increasing chain of radical polynomial ideals in C[y] asso¬ 
ciated to a Pfaffian system as in (5) (or its dual counterpart: a descending 
chain of algebraic varieties in C m ) which allows us to give a geometrical 
criterion for the solvability of the system. 

Definition 7 Let (3 p ) pe No die sequence of radical polynomial ideals in 
C[y] defined recursively as follows: 


10 



















3o := \J (g), where the equations g = 0 define the algebraic constraint 

o/S. 


• 3 1 := y £ + 3 0 , where J C C[y] is the ideal generated by the Frobenius 
conditions as in Notation 2. 



• For p> 1, 3 p+1 := y%. 

Here, for p > 0, 3 P C C[y] denotes the prolongation ideal introduced in 
Definition f. 

By duality, we define for each p > 0 the set Y p C C” as the set of zeros of 
the ideal 3 P in the affine space C n . 

We obtain an ascending chain of polynomial ideals 3o Q • • • C 3 P C 
3 P+ 1 C ••• and the corresponding descending chain of algebraic varieties 


V 0 2 ■ ■ ■ 2 V p D V p+ i D ■ ■ ■ . 


Each step 3 P 3 p+ \ or Y p V p+ i is called a prolongation. 

From Proposition 5, both chains (J p ) pe N 0 and (V p ) pe Ho depend only on 
the Pfaffian system £ and, by Noetherianity, there exists a minimum integer 
Poo where the prolongations become stationary. Moreover, from the defini¬ 
tions of the chains we observe that if an integer p satisfies 3 P = 3 p +i (resp. 
Y p = Vp+i) therefore the chain becomes stationary and hence p > poo■ In 
other words, p^ is the smallest integer p > 0 such that 3 P = 3 p +i (resp. 
V p = V p+ i) holds. 

This section is devoted to state integrability criteria for autonomous 
differential-algebraic Pfaffian systems in terms of the prolongation varieties 
Y p introduced above. 

As a consequence of Lemma 6 we deduce: 

Lemma 8 Let £ be an autonomous differential-algebraic Pfaffian system: 



where f and g are finite sets of polynomials in C [y]. Let p € No be such 
that Y p and V p+ i are the same nonempty variety in a neighborhood of a 
common regular point Q E C n . Then £ is integrable. Moreover, there exists 
a solution of £ passing through the point Q and contained in V p +i. 

Proof. Let g p C C[y] be a finite system of generators of 3 P , the defining 
ideal of V p . Consider the autonomous differential-algebraic Pfaffian system 
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Note that the system X p verifies the assumptions of Lemma 6 with X = J p 
and J = J p +i- Then X p is integrable and there exists a solution of X p 
passing through Q and contained in Vp+i. The lemma follows from the fact 
that (g) C (g p ). ■ 

Corollary 9 If the system X is not integrable then for all integer p such 
that Y p 0, the inequality dimVp + i < dimVp holds. 

Proof. If the relation dimVp = dimV p +i holds for some integer p > 0 
with Y p 0, since V p+ i C Y p , the varieties Y p and Vp+i share a common 
(nonempty) irreducible component C of maximal dimension. Then, any reg¬ 
ular point Q G C verifies the hypotheses of Lemma 6 and, therefore, X is 
integrable. ■ 

The following proposition, which is an immediate consequence of Corol¬ 
lary 9, gives a first criterion for the integrability of X: 

Proposition 10 Let X be an autonomous differential-algebraic Pfaffian sys¬ 
tem. Then, X is integrable if and only ifY Poo 0. 

Proof. First observe that the condition Y Poo 0 is clearly necessary because 
if 7 is a solution of X, its image is contained in any variety Y p , in particular, 
in V Poo and then, Y Poo 0. Reciprocally, if V Poo 0 but X is not integrable, 
Corollary 9 implies that the decreasing chain is not stationary at level p 0c , 
which contradicts the definition of the integer p^. The proposition follows. 


Remark 11 In the language of Proposition 10, the classical Frobenius The¬ 
orem can be stated as follows: assume that no algebraic constraints appear 
in the system X, then X is completely integrable if and only if Y Poo = C n 
(or equivalently 3 Poo = 0). Obviously in this case p^ = 0. 

Example 12 Consider the following autonomous Pfaffian system for m = 
n = 2 with no algebraic constraints: 


' dyi 


dy 1 

2 

dxi 

= yi 

dx 2 

= V\ 

dy 2 

= yw2 + 1 

dfm 

2 

, dx\ 

dx2 

= VI 


In this example, Jo = 0. The Frobenius conditions are y\ and y\(y 2 +Vi — 2); 
therefore, the polynomial ideal $ C C[yi,y2] is the principal ideal generated 
by y\ and hence its radical is generated by y\. In order to obtain the 
ideal J 2 , we compute the polynomials Dffyi) and ^(yi) obtaining y\ and 
y\, respectively. Then J 2 = y (yi, yf) = (y 1 ) = 3i and so, p 0Q = 1. Since 
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(t/i) / C[y 1 , 2 / 2 ] the system £ is integrable (but not completely integrablej. 
Moreover, any solution 7 := ( 71 , 72 ) must satisfy 71 ( 37 , 37 ) = 0, because its 
image must be included in the line {y\ = 0}. Hence, 72 does not depend on 

d r (2 

37 and satisfies —— = 1 and therefore 72 ( 37 ) = 37 +A for a suitable constant 

OX 1 

A G C. In other words, the solutions ofT are 7 ( 37 , 3 : 2 ) = (0 ,37 + A). 

Proposition 10 gives a conceptually simple criterion to decide the inte- 
grability of a Pfaffian system £: it suffices to compute and to check if 
the algebraic variety Y Poo is empty or not. Even if the integer p^ exists by 
Noetherianity, a priori it could be too big. However we will show that p^ 
is in fact bounded by the dimension of the ambient space. 

Proposition 13 Let £ be an autonomous differential-algebraic Pfaffian sys¬ 
tem. If, for some po G No, C C V Po is an irreducible component which is also 
included in V po+ i> then C is an irreducible component ofY Pao . Moreover, 
Y Poo = Y n+ i and then, Poo < n + 1. 

Proof. Due to Lemma 8, we can guarantee that for any regular point 
Q G C passes locally the image of a solution of £. Since every solution of £ 
is contained in all the varieties Y p , we conclude in particular that the closure 
of Reg(C) is contained in V Poo . On the other hand, the regular points form 
a dense subset of C, so C is an irreducible algebraic set contained in V Poo . 
Moreover, it is also an irreducible component of V Poo because by assumption 
it is an irreducible component of Y p and V Poo C Y p . 

If Y p is empty for some index p, then Y Poo is empty too. Therefore, 
without loss of generality we may suppose V n+ i / 0. Let C n+ \ be an 
irreducible component of V n +i. There exists a decreasing chain of irreducible 
varieties C p , p = 0,. .. ,n + 1, where each C p is an irreducible component 
of Y p . If the sequence (C p ) p is strictly decreasing, we have dim(C ri _)_i) < 
dim(C' n ) < • • • < dirn(Ci) < dim(Co) < n and then C n+ 1 = 0, leading to a 
contradiction. Hence the sequence (C p ) p stabilizes in C n+ \ for some index 
p < n. Thus, C n+ 1 is an irreducible component of Y Poo . Since this holds 
for any irreducible component of V n+ i, we deduce that V n+ i C V Poo . The 
other inclusion is always true; therefore, we have V Poo = V n +i- ■ 

Remark 14 The variety M := Reg(V Poo ) C C n is the integral submanifold 
associated to the system £ in the usual sense: for each point Q G M, its 
tangent space TqM is spanned by the vector fields determined by the differ¬ 
ential system £, and the dimension of M is maximal with this property. For 
instance, in the previous example we have Y Poo = {0} x C and the solutions 
are 7 ( 37 , x 2 ) = (0 ,37 + A). 

Hence, the invariant g := dim(V Poo ) is a measure of the integrability of 
the system £; it agrees with the maximal dimension of integral submanifolds 
of £. The extreme cases g = n, and g = — 1 correspond to the complete 
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integrability and the inconsistency o/£, respectively. In the previous example 
we have g = 1. 

We summarize Propositions 10 and 13 in the following criterion concern¬ 
ing the integrability of £: 

Theorem 15 (Criterion for the integrability of a Pfaffian system) 

Let £ be an autonomous differential-algebraic Pfaffian system. Then 

£ is integrable ■<=>■ V Poo 7 ^ 0 V n+ i / 0. ■ 


For a non-autonomous differential-algebraic Pfaffian system £, Theorem 
1 stated in the Introduction is proved by considering the associated au¬ 
tonomous system £ au t defined as in (4), and making all the constructions in 
this section from £ aut . 

3 Quantitative Aspects 

3.1 Some tools from effective commutative algebra and al¬ 
gebraic geometry 

Throughout this section we will apply some results from effective commu¬ 
tative algebra and algebraic geometry. We recall them here in the precise 
formulations we will use. 

One of the results we will apply is an effective version of the strong 
Hilbert’s Nullstellensatz (see for instance [12, Theorem 1.3]): 

Proposition 16 Let f\, . . . , f s € C[zi,...,z n \ be polynomials of degrees 
bounded by d, and let I = (f\,..., f s ) C C[^i,..., z n j. Then (\/7) dn C /. 

In order to obtain upper bounds for the number and degrees of generators 
of the radical of a polynomial ideal, we will apply the following estimates, 
which follow from the algorithm presented in [17, Section 4] (see also [15], 
[14]) and estimates for the number and degrees of polynomials involved in 
Grobner basis computations (see, for instance, [5], [19] and [7]): 

Proposition 17 Let I = (/ 1 ,..., f s ) C C[zi,..., z n ] be an ideal generated 
by s polynomials of degrees at most d that define an algebraic variety of 
dimension r and let v = max{l,r}. Then, the radical ideal \fl can be 
generated by ( sd) 2 ° ( ' 1 ' n) polynomials of degrees at most ( sd) 2 ° < ' vn) . 
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3.2 An effective decision method 


The integrability criteria presented in Section 2.5 enable the application of 
tools from effective algebraic geometry in order to derive a decision method 
for the PfafRan systems under consideration. 

As stated in Theorem 15, we have that for an autonomous differential 
system X in n differential unknowns, X is integrable if and only if V n +i ^ 0. 
We start by estimating the number and degrees of polynomials generating 
the intermediate ideals 3 p and the complexity of computing them. 

Lemma 18 With our previous notation, let v := max{l, dim(Vo)}, o := 
max{l,s} and d := max{deg(f),deg(g)}. There exists a universal con¬ 
stant c > 0 such that for each 0 < p < p^, the ideal 3 P can be generated 
by a family of polynomials g p whose number and degrees are bounded by 
(nmad) 2C v . These polynomials can be computed algorithmically within 
complexity ( nmad) 2 ° < ' n ^. 

Proof. For p = 0, we have that go is a set of generators of \J (g); then, by 
Proposition 17, they can be chosen to be at most (ad) 2 ' 0 ' " polynomials of 
degrees bounded by (ad) 2C °'' n , where Co is a universal constant. 

The ideal 3\ is the radical of the polynomial ideal generated by go, 
Dj(g) for every g in go and 1 < j < rri, and the Frobenius conditions 
Dj(fik) — Dk(fij) for all indices i,j, k with i = 1,..., n, and j,k = 1 ,m. 
For each polynomial g in go, we add m polynomials Dj(g), for 1 < j < m, 
and, since deg(f) < d, we have that deg (Dj(g)) < deg (g) + d. The Frobenius 
conditions are given by n(™) polynomials of degrees bounded by 2d. We 
may assume that the constant Co is chosen so that the number and degrees of 
all these polynomials is bounded by (nmad) 2Cavn . Then, by Proposition 17, 
the radical ideal 3i can be generated by a set of polynomials whose number 
and degrees are bounded by ( nmvrd) 2 '' co '' n+1 < ( nmad) 2,20 ’'" by taking, for 
instance, c > Co + 1. 

Assume that, for p > 1, there is a system g p of polynomials that generate 
3 p whose number and degrees are bounded by (nmcrd) 2C(P+1)l "". Recalling 

that 3 p+ i = \Jjp = sj(g p ; Dj(g p ), 1 < j < m), the bounds from Proposition 
17 imply that can be generated by a set of polynomials whose number 
and degrees are bounded by 

\ 2 c 0 l/n 

/ 2 c (p+ 1 ) I/71 ' \ 2 \ - / .2 c (p+ 1 ) l ' n + c o I ' n -+2 , 2 c (v +2 '> vn 

\[m+l){\nmad) ) ) < (nmcrd) < (nmad) 

for a universal constant c (it suffices to take c > Co + 2). 

The complexity bound follows from the complexity of the computation 
of the radical of a polynomial ideal stated in [17, Section 4]. ■ 

As a consequence of the previous lemma, we deduce that we can obtain a 
Grobner basis of the defining ideal of V n +i within complexity ( nmcrd) 20( ' n \ 
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which enables us to decide immediately whether this variety is empty or not. 
We conclude: 


Theorem 19 Let £ be an autonomous differential-algebraic Pfaffian sys¬ 
tem: 

( dy 


£ = 


dx = f(y) ' 

g(y) = o 


where y = yi,...,y n and x = xi,...,x m , f = (/ij)i<i<n,i<i< m and g = 
gi ,..., g s polynomials in C[y] with deg(f), deg(g) < d. There is a determin¬ 
istic algorithm that decides whether the system £ is integrable or not within 
complexity ( nmad) 2 ° < ' n \ where a = max{l,s}. 


3.3 An effective Differential Nullstellensatz for differential- 
algebraic Pfaffian systems 

We write C{y} for the polynomial ring in the infinitely many variables yi >a 
where a £ N™. This ring has naturally m many independent derivations 

for h £ C{y}, 

L/«t| l/ihf 


dh 

dxj 


dh 


^ dyi 


Vi,a+ej i 


where ej denotes the j-th vector of the canonical basis in N™. We identify 
the variable y t with y,^ for 1 <i < n. 

Consider an autonomous differential-algebraic Pfaffian system 

e = { i = f(y) ’ 

l g(y) = 0 

where f, g are polynomials in C[y] of degrees bounded by an integer d. This 
sytem induces a differential ideal in C{y} (that is, an ideal of C{y} which is 

d 

closed under the derivations ——) simply by taking the ideal generated by 

dxj 

the polynomials y ifii - fy, g k (with i = 1 ,..., n, j = 1,..., m, k = 1,..., s) 

dy 

and all their derivatives. We denote this ideal by [—- f, g]. 

By the differential Nullstellensatz (see [13]), £ has no solution if and only 

dy 

if 1 lies in the differential ideal [—-f, g] C C{y}. The aim of this section 

dx 

is to prove an upper bound for the number of derivations needed to obtain 
a representation of 1 as an element of this differential ideal assuming that 
£ is inconsistent. A previous bound on this subject for general differential 
algebraic systems is given in [8]. 


16 




For each a = (aq,..., a m ) G Nq\ let |ct| = a\ + ... + a m and for each 

k > 0, let y[ fc l be the set of all variables yi a with |a| < k. For h G C{y}, we 

d^h d^h 

denote —-= ^ ^ n and 


3x° 


dx ^ 1 .. .dx 


h [k] = 


d^h 

dx. a 


lal < k 


We consider the algebraic ideals 3 p C C[y] (0 < p < p^) introduced 
in Notation 7. For p = 0, ... ,Poo, let g p C C[y] be a system of generators 
of the ideal 3 p . We also denote g_i := g. Since (g) C 3 P for every p, the 
differential-algebraic Pfaffian system 


={§= f(y) - 
l g P (y) = o 


dy 

has no solution and so, by the differential Nullstellensatz, 1 G [—- f,g p ]. 

Thus, there exists a non negative integer k (depending on p) such that 
1 € ((^" — f)^,gp^) C[yt fc+1 l]. Moreover, by Proposition 10, we have 
that {g Poo = 0} = V Poo = 0, and so, 1 € (g Poo ) C C[y]. We define 

k p = min |k G N 0 : 1 € - f) [fcl , g^' 1 ) }. (10) 


In particular, 

k- 1 = min j k G N 0 : 1 G ((|^ - f) [fe] , g W ) } 

is the order of differentiation of the input equations we want to bound. 

Note that, since (g p ) = 3 P C 3 p +i = (g p +i) for every p , the sequence 
k p is decreasing. In addition, as 1 G (g Poo ), we have that k Poo = 0. We will 
obtain an upper bound for k_\ by recursively computing upper bounds for 
k p for p = p ^,..., 0. In order to do this, we will use the following auxiliary 
sequence of non-negative integers defined for p = 0 ..., p Q0 : 


e 0 : = min{e G N : C (g)}, 

£\ := min{e G N : Jf C 3o + 5}, 

£ p := min{e G N : J p C 3 p _i}, for 2 < p < p 00 . 

By the definition of the ideals 3 p and 3 p , it follows that £o is the Noether 
exponent of the ideal (g), E\ is the Noether exponent of the ideal Jo + 3 
and, for p > 2, e p is the Noether exponent of 3 p _i. Taking into account that 
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- f) [1] ) and > for P > 1, ^p-i = (gp-i;-Dj(g p -i), 1 < j < m) and 


^j(gp-i) = 


_ ggp-i 
dxj 


mod 


/% _ 


f), it follows that 


V C (g_ 


and U p p C 


V<9x 

i). •’t c (<^ _f > |I1 ’ g » 1 )’ 

f, gj'M for p> 2. 


/fly _ f „w 


VOx 


( 11 ) 


The following is the key result that enables us to prove the required 
bound for the integers k p : 

Lemma 20 Let 0 < p < p^ and g G J p . Then, for every k € No, all the 

Q\ a \ g 

partial derivatives o , for a € (No) m with |a| < k, lie in the polynomial 


ideal ^ ((|| - f)[ 1 + e p fe l,g| ) 1 J t " 1 £pfc] 


Proof. Consider the usual graded lexicographic order in ((No) m ,-<): for 
a,j3 € (No)” 1 , a -< (3 if and only if either |a| < |/3| or |a| = |/3| and there 
exists j, 1 < j < m, such that ak = f3k for every 1 < k < j and oij < f3j. 
We proceed inductively following this order. 

By the inclusions in (11), we have that g £p G ~ O^igp-i) and so, 
the statement holds for k = 0. 

Claim: Given e G N, for every a G (No) m , there exist c € N and a differen¬ 
tial polynomial G G C{y} in the ideal generated by the partial derivatives 

(Wg\ 


Ox- 


such that 


/ p<a. 


0 £ l Q l g 
0x e ' c 


= C ' 


0x Q 


+ G. 


( 12 ) 


Proof of the claim. According to the Leibniz formula, we have that 


0 £ l“l( 

dx £C 


= E« 


' d^g 
Ox^i 


’ d\ v *\ g 
dx. Ue 


where c v > 0 and v runs over all matrices in (No) exm , with rows v\,... ,v e G 
(No) m , such that the sums of their columns are ea \,..., ea m respectively. 
Let us analyze each term of the above sum. If, for some i = 1,... ,e, we 


have that |zzj| < |a|, the term lies in the ideal generated by 
it will contribute to the differential polynomial G. 


Ox/ 3 


and 


' [t-<CL 
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Assume now that \vi\ > |a| for every i. Then, \vi\ = a for every i, since 
l^il + • • • + \v e \ = e|a|. If t'i -< a for some i, the term will contribute to 
the polynomial G. Then, it remains to consider the case where no zv, is 
smaller than a in the order -<. We claim that in this case ig = a for every i. 
Otherwise, let io be such that Vi 0 is the smallest among all Vi different from 
a, and jo = min {7 : a,j < 1 %j}. By the minimality of Vi 0 and jo, all the 
entries Vij 0 of the joth column of the matrix v are greater than or equal to aj 0 
and for at least one index, io, the strict inequality holds; therefore, the sum 
of the joth column of v is greater than £otj 0 , which leads to a contradiction. 
Therefore, in the only term not containing a factor with z/j -< a we have 

/ 9l Q l g X £ 


Vi = a for every i and so, the term equals c 
constant c. 


3x° 


for a combinatorial 


Since g £ ? € ((|^ - f) [1] , 


gi-i 


it follows that 


Qe v \a\ 


g £p 


dx. £ v c 


((-- 

Qe p \a\ 


g £p 


f\[i+o P fc] [i+opA,-]\ £ (Jr ever y a g w ith |ct| < k\ in particular, 

' p / ' ax £ p' a 

lies in the radical of this polynomial ideal. By induction in ((No) m , A), this 

property along with formula (12) of the Claim imply that 


gH 


g 


3x° 




[1+Epfc] [1+Epfc] 

’ Op—1 


for every a € (No) m with |a| < k. ■ 

Corollary 21 LetO<p<poo. For every k € No, if 1 € ~ f)^, , 

then 1 € — f)[ 1 + e P fe ] j gj ) 1 J t " 1 Epfc ^. In particular, k p -\ < 1 + e p k p . 

Proof. Since e„ > 1, it follows that — f)^ C — f)[ 1+e * fc ]. In 

ox Ox 

addition, applying Lemma 20 to the polynomials in g p , we deduce that 
(gf 1 ) C j ((|| Then, 

((% - o'* 1 .*? 1 ) c ]/{(%-f) ,1+e - tl 4F kl )- 

which implies the first assertion of the Corollary. 

The second assertion follows from the definition of the integers k p in 
equation (10). ■ 


M-i 

Corollary 22 Let p := min{0 < p < Poo : k p = 0}. Then, k-\ < (p + 1) f~[ £ v 

p =0 
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Proof. By applying recursively the inequality k p -\ < 1 + e p k p from the pre- 

vious corollary, it follows easily that k p -j < j n £i for j = l,...,n,H + 

i=p-j +1 


1 . ■ 


Now, to prove our main result, it suffices only to bound the Noether 
exponents e p for p = 0,— 1, which will be easily done from Proposition 
16. 

Theorem 23 Let x = aq,.. ., x m independent variables and y = yi, ■ ■ ■ ,y m 
differential variables, and let f = ( fij)i<i<n,i<j<m and g = gi,...,g s be 
polynomials in C[y] of degrees bounded by d. Let V C C n be the variety 
defined as the set of zeros of the ideal (g) and u := max{l, dim(V)}. Then, 



C z 

where K < ((s + nm 2 )d) 2 " " for a universal constant C > 0. 

Proof. Following Corollaries 22 and 9, notice first that g < p^ < v + 1. 
On the other hand, by Proposition 16 and Lemma 18, it follows that the 
Noether exponents e p can be bounded as follows: 

£o < d n and e p < ( nmad) n2 (p+ ) , 1 < p < poo. 

By Corollary 22, we have that, for a suitable constant c, 

I\ = i < (n + 2) JJ(n7TUTd) n2c(P+1>I " n < 

p =0 

<{n + 2){nmad) < (nmad) 2 

for a suitable constant C. ■ 
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